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I determine the scaling behavior of the free energy barriers encountered by 
a flux line in moving through a three-dimensional random potential. A com- 



^ . bination of numerical simulations and analytic arguments suggest that these 

o , 

. barriers scale with the length of the line in the same way as the fluctuation 

. in the free energy. 
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Magnetic flux lines (FL) in high-Tc superconductors are one of the simplest examples of 
glassy systems In thermal equilibrium, a FL is pinned by defects (oxygen impurities, 

grain boundaries, etc) in the superconductor which lower its energy . This effect is limited 
by the line tension which opposes the bending of the line. The resulting free energy landscape 
for the FL is rather complicated and has many local minima, i.e. metastable states [Q. When 
an electric current flows through the system, the FL feels a Lorentz force perpendicular to 
its orientation and to the current direction. As long as the current is not strong enough to 
overcome the pinning forces, the line moves by thermally activated jumps of line segments 
between metastable conflgurations 1^-0]. The length of these line segments is estimated by 
the condition that the free energy barrier for a jump should be of the same order as the 
gain in free energy due to that jump. These dynamics are believed to be the reason for the 
nonlinear volt age- current characteristics found in experiments 0. 

Since energy barriers play such an important role in the dynamics of FLs, it is essential 
to know their properties. The scale of these barriers should grow with observation size L 
like a power law . Usually, it is assumed that the energy scale in the system is set by the 
fluctuations in free energy which increase as L^, and that therefore ip = 6 1^,0]. However, 
it is also quite possible that the heights of the ridges in the random energy landscape scale 
differently from those of the valleys that they separate, with ip > 9. Yet another scenario is 
that transport occurs mainly along a percolating channel of exceptionally low energy valleys 
with ip < 6. A flrst attempt to clarify this situation was taken in |p, where ip = 6 was 
established for a FL moving in 2 dimensions. This demonstration relied on exact results 
for minimal energies in 1 + 1 dimensions, and on the fact that the endpoint of a FL in 2 
dimensions has to move through all points which lie between its initial and flnal positions. 
It is of great importance to discuss also a FL in 3 dimensions, which is more physically 
relevant. In contrast to a two-dimensional system, a FL which moves in three dimensions 
can avoid regions in space which are energetically unfavorable for one of its segments, and one 
might therefore speculate that ip < 6 in three dimensions. In this paper, we flrst determine 
numerically a lower bound for the barrier energy which scales in the same way as the energy 
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fluctuations, thus ruling out ip < 6. Further numerical results predict that an upper bound 
scales in the same way, ruling out ip > 0, and thus leading to ip = 9. 

We describe the FL as a directed path in a random medium |Ty] . The path is discretized 



to lie on the bonds of a cubic lattice, starting at the origin and directed along its (1,1,1) 
diagonal. Each segment of the line can proceed in the positive direction along one of the 
three axes, leading to a total of 3* configurations after t steps, with endpoints lying in the 
plane which is spanned by the points (t,0,0), (0,t,0), and (0,0, t). A given configuration 
of the FL is labelled by vectors {x(t)} for r = 0, 1, ■ ■ ■ , t, giving the transverse coordinates 
of the FL at each step. The points {x{t)} lie on the vertices of a triangular lattice. For a 
given value of r, they lie on one of three alternating sublattices. 

To each bond on the cubic lattice is assigned a (quenched) random energy equally dis- 
tributed between and 1. The energy of each configuration is the sum of all random bond 
energies on the line. For each endpoint {t,x), there is a configuration of minimal energy 
Emin{x\t) which cau be obtained numerically in a time of order by a transfer matrix 
algorithm |T^. The fluctuation in minimal energy is known to scale as with 6 ~ 0.24, 
and the transverse fluctuation of the coordinates of minimal paths is known to scale as t^. 



with ^ 0.62 [O. The endpoints of the minimal paths with the lowest energy lie within a 
distance oc if:'' of the origin. Fig. [l| shows the minimal energies of paths of length t = 288 to 
endpoints x with |x| < 0(t''). The highest energy in this region is represented in white, the 
smallest energy in black. The minimal energies are correlated over a distance of the order 
of t^. The distribution of minimal energies is close to a Gaussian and is shown in Fig. |^. 



Similar to a 2-dimensional system [|T^, this distribution seems to have a third cumulant 
since it is not completely symmetric. 

We next examine the energy barrier that has to be overcome when the line is moved from 
an initial minimal energy configuration between (0, 0) and {t, Xj) to a final one between (0, 0) 
and {t,Xf), with \xf^i\ < t^. The only elementary move allowed is flipping a kink along the 
line. Thus the point (r, x) can be shifted to (r, x ± Cj), where ±ej are the six vectors which 
connect a vertex in the triangular lattice to its nearest neighbors within the same sublattice. 



Each route from the initial to the final configuration is obtained by a sequence of such 
elementary moves. For each sequence, there is an intermediate configuration of maximum 
energy, and a barrier which is the difference between this maximum and the initial energy. 
In a system at temperature T, the probability that the FL chooses a sequence which crosses 
a barrier of height Eb is proportional to ex-pl—Es/T), multiplied by the number of such 
sequences. We assume that, as is the case for the equilibrium FL, the "entropic" factor of 
the number of paths does not modify scaling behavior. Thus at sufficiently low temperatures 
the FL chooses the optimal sequence which has to overcome the least energy, and the overall 
barrier is the minimum of barrier energies of all sequences. 

Since the number of elementary moves scales roughly as the volume of a cone which 
contains the initial and final fines, the number of possible sequences grows as t^'^. This 
exponential growth makes it practically impossible to find the barrier by examining all 
possible sequences, hampering a systematic examination of barrier energies. Rather than 
finding the true barrier energy, we proceed by placing lower and upper bounds on it. 

A lower bound to the barrier energy is obtained in the following way: While the line moves 
from its initial to its final configuration, the transverse coordinates of its endpoint move 
between nearest-neighbor positions on one of the above mentioned triangular sublattices. 
When the endpoint is at a position (x), the energy of the line is at least as large as the 
minimal energy Emin{x\t). The maximum of all these minimal energies along the trajectory 
of the endpoint, minus the energy of the initial configuration, certainly bounds the barrier 
energy from below. Since wc do not know the actual trajectory of the endpoint, we have to 
look for the trajectory with the smallest maximal energy. Only in this case we can be sure 
that we have indeed found a lower bound. This situation is fundamentally different from a 
2-dimcnsional system, where there is only one possible trajectory for the endpoint. 

Provided that the minimal energies E^in{^\t) are known, this lower bound is determined 
in polynomial time by using a transfer-matrix method: We start by assigning to the initial 
point Xi a lower bound energy 0, and to all other sites x on the same sublattice an energy 
t which is certainly larger than the lower bound resulting from the algorithm after many 
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iterations. At each step the energy of all sites x except the initial site is updated according 
to the following rule: Look for the minimum of the energies of the 6 neighbors 5 ± e,. 
If this is smaller than the energy at x, replace the energy at x by this minimum or by 
Emin{x\t) — Emin{xi\t), whichcver is larger. After a sufficiently large number of iterations, 
which is of the order of the size of the area of interest (which scales as t^''), all possible 
trajectories to endpoints within this area have been probed, and the energies do not change 
any more. The energy at site Xf is then identified as the lower bound. Fig. ^ shows the 
lower bound to the energy barrier for a line with the endpoint moving from the origin to 
sites within a distance of the order of t^, for different values of t and averaged over 500 
realizations of randomness. The distance \xf — Xi\ has been scaled by t"'', and the energy 
by t^^. With this scaling, all the curves should collapse, leading to the following scaling 
behavior for the lower bound, 

{E^{t, \xf - m = t'Ulxf - x,\/t^). (1) 

The function f{y) is proportional to y^^'^ for small y. For the simulated system sizes, however, 
this asymptotic scaling cannot yet be clearly seen. For y > 1, the scaling form in eq. (|ip 
breaks down since the minimal energy is then a function of the angle {\x\/t). We conclude 
that the lower bound to the barrier scales in the same way as the fluctuations in minimal 
energy, and consequently the energy barrier increases at least as t^, leading to ip > 6. The 
distribution P{EJ) of the lower bound energy for a given distance \x\ oct^ is shown in Fig. ^. 
It appears to be half-Gaussian with width oc t'^. 

The result ip > 6 is not surprising if one realizes that an even simpler lower bound is 
given by max{Emin{x f\t) — Emin{xi\t), 0), which evidently scales as since the distribution 
function of minimal energies decays exponentially fast, i.e. has no power-law tails (see 
Fig. 0). To make sure that the scaling of the lower bound found above is not dominated 
by the neighborhood of flnal conflgurations with particularly high energies, I repeated the 
above simulations by allowing only endpoints with minimal energies smaller than the initial 
energy. This corresponds to the situation that the endpoint of the line does not move to an 
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arbitrary position but to a position which is energetically favorable. The result is shown in 
Fig. ^ and has the scaling form 

{E4t, \xf - x,\)) = t'f4\xf - x.\/t^). (2) 

As in the previous case, the asymptotic scaling /_ (y) oc y^^'^ for small y cannot yet be clearly 
seen. The energy distribution of the lower bound is again a half-Gaussian of width oc if:'' and 
looks similar to Fig. |. 

The same scaling behavior is also found when instead of the optimal trajectory for the 
endpoint the shortest trajectory (a straight line) is chosen. In this case, the mean of the 
barrier energy Eq has the scaling form 

{Eo{t, \xf - Xi\)) = t'foilxf - x,\/t<) (3) 

(see Fig. |), again with a half-Gaussian distribution of width oc t^. This, of course, does not 
represent a lower bound to the true barrier, but it will be important for the determination 
of an upper bound below, and is therefore included here. 

The result E^ oc (eq. (^) can be explained from the exponential tails of the distribu- 
tion of minimal energies: If we asume that the endpoint of the line moves only in valleys of 
particularly low energy, we can successively remove all sites with the largest minimal energy 
from the set of possible endpoints, until the connectivity over the distance breaks down. 
The remaining endpoints form percolation clusters, and their density is given by the corre- 
sponding percolation threshold (This is analogous to random resistor networks describing 
the hopping resistivity for strongly localized electrons. The resistance of the whole sample is 
governed by the critical resistor that makes the network percolate Since the occupied 

sites are correlated over the distances considered, the value for the threshold is different 
from the site percolation threshold of 0.5 in an infinite triangular lattice with no correlation 
between occupied sites. But for the present purpose, it is sufficient to know that this thresh- 
old is finite and that therefore a finite percentage of all sites are below threshold. Since the 
distribution of minimal energies dacays fast, its tail cannot contain a finite percentage of all 
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sites. We conclude that the threshold is within a distance of t from the peak, and therefore 
that the energy fluctuation on the percolation cluster and consequently the lower bound for 
the barrier are oc . 

We now proceed to construct an upper bound to the energy barrier. To this purpose, 
we specify a sequence of elementary moves which take the line from its initial to its final 
configuration. Since we cannot be sure that this sequence is the optimal one, we know only 
that the barrier associated with this specific sequence is an upper bound to the true barrier. 
The algorithm for the motion of the line is inspired by the one presented in and is as 
follows: First, one choses a sequence of endpoints connecting the initial to the final endpoint 
which is as short as possible. Then, one draws all the minimal paths leading to these 
endpoints. It is certainly advantageous to keep the intermediate paths as close to minimal 
configurations as possible and therefore to require that the line passes successively through 
all these intermediate minimal configurations. Usually, minimal configurations {xi{t)} and 
{^2(''")} with neighboring endpoints have large parts in common and separate only during 
the last few steps. They enclose a small loop with a size of the order of 1. But sometimes, 
both paths already separate during the first few steps and form a large loop of the lateral 
size of the order of t^. We have to give a prescription for how the line moves over a loop. If 
the two minimal paths have nowhere a distance larger than 1 (measured in units of |ej|), we 
can choose a sequence of elementary moves such that at most two bonds of the line are not 
on one or the other minimal path, leading to a barrier of order 1 between the two. If the 
distance is larger than 1, we proceed as follows: Let tq be the last point which both lines 
have in common, i.e. Xi(ro) = X2(to) and Xi{tq + 1) 7^ 3^2(^0 + !)• We then consider the 
midway points (tq + {t — to)/2,x) which connect both lines in the middle of the loop via 
the shortest possible trajectory (if there are several possibilities, we choose one at random). 
For each of these points, we find two minimal segments of length {t — Tq)/2 connecting on 
one side to (To,a;i(ro)) and on the other to either {t,xi(t)) or (t,X2(t)). There are usually 
several possibilities in making these connections, but all of them lead to the same result. 
The two segments form an almost minimal path of length t — tq, constrained to go through 

7 



the point (tq + (t — to)/2,x). We next move the hne {^i(t)} with tq < t < t stepwise 
through this sequence of almost minimal paths. At each step we first attempt to move 
the upper segment and then the lower one. The prescription for moving these segments 
of length (t — To)/2 is exactly the same as for paths of length t: If the distance between 
two consecutive configurations is larger than 1 for some r, we consider the points in the 
middle of the loop formed by the two, and construct minimal paths of half the loop length 
connecting them to the initial and final loop points. Next we attempt to move segments 
of the length of the loop by repeatedly moving the upper and lower line portion. In some 
cases, it is necessary to proceed with this construction until the cutoff scale (1) is reached. 
Thus, at each intermediate configuration the line is composed of segments of minimal paths 
of different length, the smallest segment having length 1 in the worst case. 

We now estimate the barrier energy resulting from the above construction. In principle, 
this can be done by programming the algorithm and determining the result numerically. 
Such a program would be more complicated than in 2 dimensions, and it would not be able 
to simulate large systems, since a big portion of the bond energies on the three-dimensional 
lattice need to be stored. It is therefore uncertain if the asymptotic scaling behavior of the 
upper bound may be found this way. Instead, we resort to analytic considerations: Since 
the line is always composed of segments of minimal paths, whose scaling properties are 
known, we have enough information to give an upper bound to the barrier energy. For small 
distances \xf — Xi\, the initial and final path usually differ only in the last few steps, and 
therefore the upper bound to the barrier energy (and also the barrier energy itself) does 
not depend on t, but only on \xf — Xi\. We are interested in the barrier which has to be 
overcome when the line moves over a distance of the order of t^. Let us add successively 
the contributions to the upper bound which result from the different steps in the algorithm: 
The first step consists in finding a sequence of minimal paths with endpoints lying on the 
shortest trajectory from Xi{t) to Xf{t). The energy difference between the minimal path with 
highest energy and the initial minimal path is EQ{t,t'^) ~ , which is the first contribution 
to the upper bound. When moving from one minimal path to the next, the line has to 
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overcome a loop which in the worst case has the length t and which might (again in the 
worst case) occur in combination with the minimal path with highest energy. We therefore 
have to add to the upper bound the contribution of a loop of length t and width cx f^. This 
is obtained as follows: Within the loop, the line moves through a sequence of paths which 
are composed of two pieces of minimal path of half the looplength. Both the upper and the 
lower sequence cover an energy range EQ{t/2, ~ (l/2)^£'o(t, t^). In the worst case, 

both sequences have their maximum simultaneously, giving a contribution (l/2)^~^£'o(i, i^) 
to the upper bound. All further contributions can immediately be written down because of 
the recursive definition of the algorithm: While the upper and lower segments move through 
minimal configurations within the loop, they in turn have to overcome loops which in the 
worst case have the size and so on. The sum of all these contributions, averaged over 
different realizations of randomness, is 

{E,{t, t^)) = {Eo{t, t^)) + 2 {Eoit, t^)) ((1/2)^ + (1/4)^ + . . . 

= {Eo{t,t<)) (2/(1 - (1/2)^ - 1) ^ 12.0{Eo(t,t<)). (4) 

In principle, one has to add a constant which accounts for the breakdown of the scaling form 
of the energy increase for small loops. But this constant is of the order of one and can be 
neglected with respect to the terms which increase with t^. 

There are several configurations of the path which are expected to have the energy Ec- 
They pass through loops of all sizes and are composed of one minimal segment of length 
t/2, one of length t/4, etc; ending with two smallest pieces of length 1. Each of these is 
a candidate to be the barrier path in our algorithm. To obtain the mean value for the 
maximum of their energies we need to know their number and their distribution, especially 
in the large-energy tail. The exact number of candidate barriers is not known, but we 
can be sure that it increases no faster than t^~^^ which is the order of the total number of 
intermediate configurations of the path. The energy distribution for the candidate barriers 
results from the energy distribution of their segments: Each of the segments of length has 
an approximately (half)-Gaussian energy distribution with a width of the order of t^. Since 
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the different segments are constructed through a specific recursive procedure, they might 
not be independent. If we want to be sure to estabhsh an upper bound, we have to assume 
the worst case that they are completely dependent, resulting in a variance 

var f)) ~ log2(t)(var [Eo{t, f)) + 2var (Eo{t/2, {t, 2f)) + • • •) 

~ 8.7 \n(t) var (^Eo{t, t^)) oc \n(t) t^^ . (5) 

It can be checked easily that (for large N), the maximum of N independent Gaussian 
variables of mean a and variance a^, is a gaussian of mean a + ayThiN and variance 
cr^/(21nA^). Since the candidate barriers have large segments in common, their energies are 
not independent. Assuming their independence, we overestimate again the barrier energy, 
but of course we still establish an upper bound to it. Putting all contributions together and 
taking into account the behavior for small \xf — Xi\, we finally obtain the following estimate 
for the upper bound in barrier energy: 

{E+{\xf-Xi\,t)) = {E^{\xf - I , i)) + ^2 In NYS.rE^{\xf - Xi\,t) 

c^(lnt)tV+(|f^-f,|/tC). (6) 

To conclude, 1 have shown that the energy barrier encountered by a FL moving in a Zd 
random medium has an upper and a lower bound which both increase with i^, except for 
logarithmic corrections. Prom this follows that the barrier itself increases with , confirming 
the hypothesis ip = 6. Since the argument presented in this paper is mainly based on the 
exponential tails of the minimal energy distribution, it can be expected that the result ip — 9 
holds also in higher dimensions, provided that the distribution of minimal energies still has 
exponential tails. 
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FIGURES 

FIG. 1. Minimal energies of paths of length t = 288 to endpoints x with |x| < 0{t^). White: 
High energies; Black: Low energies. 

FIG. 2. Probability distribution P{Ejnin) of minimal energies £^mm(0|144), averaged over 50000 
realizations of randomness. The solid line is a Gaussian distribution. 

FIG. 3. Scaling functions /-(y), /-(y), and /o(y) defined in eqs. (1) - (3) for t = 72 (solid), 
t = 144 (dotted), t = 288 (dashed), t = 576 (long dashed), and t = 1152 (dot-dashed), averaged 
over 500 realizations of randomness. The straight line has the slope O/C = 0.39. 

FIG. 4. Probability distribution of The parameters and symbols are the same as in Fig.3. 
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